Rashba Effect at Magnetic Metal Surfaces 
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We give experimental and theoretical evidence of the Rashba effect at the magnetic rare-earth 
metal surface Gd(0001). The Rashba effect is substantially enhanced and the Rashba parameter 
changes its sign when a metal-oxide surface layer is formed. The experimental observations are 
quantitatively described by ab initio calculations that give a detailed account of the near-surface 
charge density gradients causing the Rashba effect. Since the sign of the Rashba splitting depends 
on the magnetization direction, the findings open up new opportunities for the study of surface and 
interface magnetism. 
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PACS numbers: 71.70.Ej, 72.25.-b, 73.20.-r, 85.75.-d 

A key issue in condensed-matter research aiming at 
future spintronic devices Q is to control and manipu- 
late the electron spin in a two-dimensional electron gas 
(2DEG) of semiconductor systems without the need of 
applying an external magnetic field. Rashba had real- 
ized early on Q that this can be achieved by an electric 
field which acts as a magnetic field in the rest frame of a 
moving electron. The interaction between the spin s of a 
moving electron of momentum Kk with an electric field 
oriented along the z-axis e z is described by the Rashba 
Hamiltonian 



'Hr = afi(e 2 xk)-s. 



(1) 



The Rashba parameter an is proportional to the electric 
field and depends on the effective, material-dependent 
spin-orbit coupling (SOC) strength. In nonmagnetic sys- 
tems the Rashba effect lifts the spin-degeneracy of the en- 
ergy dispersion e(k) of an electronic state, and the energy 
difference between ej (k) and (k) is called Rashba split- 
ting Ae(k) = oiji |k|. Even though spintronic research 
currently focuses on spin-polarized electrons in semicon- 
ductors 0] , it is important to explore the Rashba effect 
in other material classes as well. 

A necessary condition for the Rashba effect to occur is 
the absence of inversion symmetry and, while in the pro- 
posed FET-type spin transistor || a gate voltage must be 
applied to break inversion symmetry of the 2DEG, this 
condition is naturally fulfilled by the structural inversion 
asymmetry (SIA) existing at any crystal surface or in- 
terface. Owing to SIA, electrons in a two-dimensional 
surface or interface state experience an effective crystal 
potential gradient perpendicular to their plane of prop- 
agation, hereby optimizing (e^ x k) in Eq. One 
should expect that the Rashba effect is a general sur- 
face and interface phenomenon, but up to now Rashba 
splittings have only been observed for surface states at 
Au(lll) 0,0] and W(110) [10. Recently relativistic 
density functional theory (DFT) calculations were able to 
reproduce the observed splitting of the Au sp-like surface 



statejlfjj and the analogy to a 2DEG has been pointed 
out Yet, it is still a challenging task to give a physi- 
cal picture of the Rashba effect from the electronic struc- 
ture point of view. 

This Letter presents the first experimental and theo- 
retical evidence of a Rashba splitting of exchange- split 
two-dimensional electron states. Using the surface state 
of ferromagnetic Gd metal as example we report on the 
novel finding of a k-dependent contribution to the bind- 
ing energy of this state that changes sign upon magne- 
tization reversal. It is further demonstrated that the 
Rashba effect is enhanced upon formation of a surface 
oxide layer. The experimental observations are quantita- 
tively described by ab initio calculations showing that the 
enhancement is caused by a substantial change of charge- 
density gradients at the interface between the surface- 
oxide layer and the bulk-like metal, which leads to an ad- 
mixture of p character to the d-derived two-dimensional 
state. 

In particular, we show that the Rashba effect in mag- 
netic systems bears interesting consequences when the 
spin-degeneracy of e(k) is already lifted by an exchange 
splitting AE ex , separating majority (f) and minority 
(J.) electrons. The dispersion of a state subject to the 
Rashba effect can then be written as e^)(k) = e(k) ± 
\ AE ex ± 5 aji |k| and the Rashba splitting Aeuf)(k) = 
£|(l)(k,M) — ejm(k, — M) = ±a_R |k| can be observed 
separately for majority and minority states by two sub- 
sequent measurements with opposite sample magnetiza- 
tion, M and M. Since in many cases the Rashba split- 
ting is non-zero only for surface and interface states, the 
observation of this quantity opens a new and powerful 
spectroscopic path to discriminate surface and interface 
magnetism from bulk magnetism. 

For the present investigation we chose the (0001) sur- 
face of ferromagnetic hep Gd metal where the exchange 
interaction is known to separate majority and minority 
branches of the d-derived surface state by several hun- 
dred meV 0| . With this substantial exchange splitting 
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FIG. 1: (color online) (a) Experimental geometry. The sam- 
ple magnetization M is perpendicular to the plane defined by 
the surface normal e z and the electron momentum k. (b,c) 
Angle-resolved PE spectra at 6° electron emission angle w.r.t. 
the surface normal (excited with p-polarized radiation), (b) 
Gd(0001) majority-spin surface state Sj (36 eV photon en- 
ergy), (c) Exchange-split p(l x l)O/Gd(0001) interface states 
S| and S| (45 eV photon energy). The peak positions change 
upon magnetization reversal owing to the Rashba effect. 



a limitation in previous studies of non-magnetic systems 
is avoided, where the Rashba splitting of a state, in order 
to be resolved, must be larger than its intrinsic (lifetime) 
width • The remanent in-plane magnetization of thin 
Gd(0001) films ^3 means that the spin-quantization axis 
s is orthogonal to the surface electric field, providing an 
advantageous geometry for observing a Rashba splitting, 
see Eq. (1). The p(l x l)O/Gd(0001) surface oxide al- 
lows furthermore to study an exchange-split pair of two- 
dimensional states 01 that are both occupied. 

Angle-resolved photoemission (PE) experiments were 
performed using display-type electron analyzers at the I- 
311 undulator beamline at MAX-Lab, Lund, and at the 
BUS beamline at BESSY, Berlin. In the experiments, 
the energy resolution was set to 30 meV and the an- 
gular resolution to 0.5°. Gd(0001) films were prepared 
in situ by metal vapor deposition on a W(110) single 
crystal mounted on a liquid- He flow cryostat. The film 
thickness was 10 nm as determined by a quartz balance. 
During evaporation the pressure in the UHV chamber 
rose from 5 x 10 -11 mbar to 5 x 10 -10 mbar. The sur- 
face monoxide layer was prepared according to the recipe 
given in Ref. |14| with 2 Langmuir oxygen exposure; it 
was controlled in situ by PE and LEED. The PE spec- 
tra were recorded from remanently magnetized samples 
in the experimental geometry shown in Fig.^a). 

For ferromagnetic Gd metal only the majority com- 
ponent St of the d-derived Gd(0001) surface state is oc- 
cupied [l2|, as is shown by the angle-resolved PE spec- 
trum in Fig. 1(b). At off- normal emission angles, the 
Sf peak position shifts significantly upon magnetization 
reversal. Derived from experimental spectra, the disper- 
sion near the center of the surface Brillouin zone along 
the F-M direction is presented in Fig. Ufa); filled and 
open data points distinguish the two branches measured 
for opposite magnetization directions (cf. Fig. 1(a)). The 
branches are shifted symmetrically with respect to k = 0, 
in agreement with the triple- vector product of Eq. (1). 
The energy difference Ae(k||) is plotted in Fig. 2(c). It 
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FIG. 2: (color online) (a) Experimental dispersion of the ma- 
jority spin Gd(0001) surface state in the T — M azimuth; solid 
(red) and open symbols correspond to opposite magnetiza- 
tion directions, see Fig. 1(a). (b) Sf (majority spin) and Sj 
(minority spin) interface states of p(l x l)O/Gd(0001) in the 
r — K azimuth, exchange split by AE ex — 450 meV at V. (c) 
Rashba splitting Ae derived from the data in (a) (indicated 
by Gd-Sf ) and (b) (indicated by 0/Gd-S T and O/Gd-Sj). 



is identified as Rashba splitting of the Gd(0001) surface 
state; it remains rather small (< 25 meV) in the kii range 
of ±0.3 A" 1 . 

The p(l x l)O/Gd(0001) surface oxide exhibits an 
exchange-split pair of surface bands Q| that are both 
occupied, shown as Sf and Sj in Fig. ^c). For nonzero 
k|| (off-normal emission), the peak positions clearly 
change into opposite directions when the magnetization 
is reversed. Their (upward) dispersions are shown in 
Fig. 2(b): they change symmetrically with respect to the 
Brillouin zone center, strongly supporting our interpre- 
tation as Rashba splitting. The energy splitting Ae(kii) 
is plotted separately for Sf and S; in Fig. 2(c); it is sub- 
stantially larger (about three times at jkyj = 0.3 A -1 ) 
than the Rashba splitting of the majority surface state 
on Gd(0001). Moreover, the Rashba parameters for the 
oxide-covered and the clean metal surface are found to 
have opposite signs, see Fig. 2(c). - We note in pass- 
ing that for p(l x l)O/Gd(0001) also the modulus of the 
slopes, | or |, is different for Sf and Sj. 

The ab initio calculations were performed using DFT 
in the local density approximation (LDA) employing the 
form of Moruzzi et al. 15]. We use the full-potential 
linearized augmented-plane-wave method in film geome- 
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FIG. 3: (color online) (a) Calculated majority-spin surface- 
state dispersions of Gd(0001) along two different high- 
symmetry directions (shown as red lines) of the magnetic 
surface Brillouin zone, given as inset. Hatched areas indi- 
cate bulk band regions. For symmetry labels, see text, (b) 
Charge density distribution (isolines on log scale) for Sj at 
r in a plane perpendicular to the surface. The linear color 
scale ranges from red (high charge density), yellow (medium), 
green (low), blue (very low) to white (negligible). 

try [IH Hrf. with SOC included self-consistently accord- 
ing to Ref. 01 ■ F° r a proper description of the 4/ elec- 
trons we applied the LDA+U method |l9j . The Gd sur- 
face was simulated by a structurally relaxed 10-layer film 
embedded in two semi-infinite vacua. A plane-wave cut- 
off of fc max = 3.8 (a.u.) -1 was used, and the irreducible 
part of the two-dimensional Brillouin-zone (BZ) was sam- 
pled at 21 special k|| points (36 k|| points for calculations 
with SOC included). 

In the calculations with SOC included, we took into ac- 
count that the in-plane magnetization lowers the p3ml 
symmetry of hep Gd(0001) to cm symmetry. One can 
still label the high symmetry points of the BZ as M and 
K, but with subscripts to distinguish between incquiva- 
lent points, see Fig. |3{a). Points that are related by an 
inversion center are primed to indicate that the band- 
structure is equivalent to the unprimed ky -points if one 
interchanges spin-up and spin-down bands. The magne- 
tization was assumed to point in the direction of the in- 
plane nearest neighbors. The size of the Rashba splitting 
decreases from the T-Mi direction to T-Ki and further 
to T-M2, and is zero in T-K2 direction. 

The theoretical results forp(l x l)O/Gd(0001) are pre- 
sented in Fig. 4. In the calculation, the O atoms were 
'adsorbed' on the fee site and its relaxed position was 
found 0.78 A above the plane of Gd surface atoms, s- 
Gd. Upon O adsorption, the s-Gd layer shows a strong 
(18%) outward relaxation, while the position of the in- 
ner layers remains almost unperturbed. Very similar re- 
sults were obtained for O adsorbed on the energetically 




FIG. 4: (color online) (a) Calculated surface-state dispersions 
of p(l x l)O/Gd(0001) along the high symmetry directions F- 
Mi (left side) and F-Ki (right side) of the surface BZ, see red 
lines in Fig. 3(a); (b) Rashba splittings, Ae, for majority and 
minority bands (triangles) and for the clean Gd surface (red 
crosses) given for comparison, (c) Charge density distribution 
for S| at F; a similar plot has been obtained for S|. The linear 
color code is like in Fig. 3. 

slightly less favorable hep site. In Fig. 4(a), the two oc- 
cupied states disperse upwards (positive effective mass) 
in good agreement with experiment, see Fig- El (The ex- 
change splitting is somewhat larger than the experimen- 
tal value, an overestimation that had been noted before 
for clean Gd(0001) |l9|.) The calculated Rashba splitting 
Ae of majority-spin and minority-spin bands is shown in 
Fig. 4(b); it is about three times larger than that of the 
Gd(0001) metal surface, and cir is found to be of oppo- 
site sign, in excellent agreement with experiment. 

With the results of the present ab initio calcula- 
tions, we can address the question why the oxide 
layer on Gd(0001) causes an enhanced Rashba splitting. 
Figs. 3(b) and 4(c) display the charge density distribu- 
tions of the two-dimensional states at T for the two sys- 
tems, Gd metal and O/Gd, in a plane perpendicular to 
the surface. In addition to the conventional isolines (log 
scale), the surface-state charge densities are also given 
on a linear scale (color) for better visibility. When com- 
paring the profiles it becomes obvious that adsorption 
of the electronegative O atom changes the charge den- 
sity distribution of the entire near-surface region. While 
the surface state of the metal surface, cf. Fig. 3(b), re- 
sides almost exclusively (to ~ 89%) in the top surface 
layer s-Gd, the two-dimensional state in p(l x 1)0/Gd 
is distributed over both the s-Gd and (s — 1)-Gd layer, 
see Fig. 4(c). Hence, given the close vicinity of O and 
s-Gd layer, this state may be conceived as interface state 
located between a p(l x 1)0/Gd surface monoxide layer 
and bulk-like Gd metal. 

At Gd(0001) the surface-state charge is located quite 



4 



symmetrically below and above the s-Gd plane, cf. 
Fig. 3(b), corresponding to a small (positive) charge den- 
sity gradient along the surface normal. The small charge 
density gradient directly indicates that there is a small 
electric field in the s-Gd layer of Gd(0001), which gives 
rise to the small Rashba splitting, shown in Figs. 2(c) 
and 4(b). At p(l x l)O/Gd(0001), by contrast, the rather 
asymmetric charge distribution of the interface state in 
the (s — 1)-Gd layer directly indicates the presence of a 
high effective electric field at this layer. The charge gra- 
dient is negative, i.e. opposite to the uncovered metal 
surface. We are led to conclude that it is this reversed 
effective electric field in the (s — 1)-Gd layer that causes 
the opposite sign of the Rashba parameter of the interface 
state inp(l x l)O/Gd(0001) as compared with Gd(OOOf). 

Furthermore, the calculations reveal that the surface 
oxide formation is accompanied by a pronounced change 
in orbital character of the two-dimensional state. In the 
s-Gd layer it changes upon oxidation from almost ex- 
clusively d z 2-like (d : p-ratio w 8 : 1, integrated over the 
muffin-tin sphere) to predominantly s-like with an admix- 
ture of other orbital components (d : p : s w 3 : 2 : 11). 
In the (s — 1)-Gd layer of p(l x 1)0/Gd the state remains 
d z 2-like yet with a substantial p z -admixture (d : p « 5 : 
1). Moreover, in order to identify the relative contribu- 
tions of the individual layers to the Rashba splitting, we 
calculated hypothetical values of the splitting with SOC 
set to zero for all other layers. Again there is a strik- 
ing difference of the two systems: at Gd(0001), the s-Gd 
layer provides by far the main contribution (~ 90%) to 
the Rashba splitting; yet in the surface oxide system it 
is the (s — 1)-Gd layer that accounts for over 70% of the 
splitting, whereas the s-Gd layer of p(l x 1)0/Gd cannot 
contribute owing to its prevailing s character (see above) . 

In the light of these results we arrive at the following 
physical picture of the Rashba effect. At the Gd(0001) 
metal surface SIA leads to a small but significant spill- 
out of the d-derived surface state, yet the charge gradient 
is small indicated by the relatively small admixture of p z - 
character (antisymmetric w.r.t. the surface plane) to the 
<i z 2-derived state. With the epitaxial p(l x 1)0/Gd sur- 
face oxide layer present, the electronegative O attracts 
charge from the s-Gd atom so that strong charge-density 
gradients arise not only in the s-Gd layer but also in the 
sub-surface layer. As a consequence, the two-dimensional 
state changes in spatial distribution (becoming an inter- 
face state) and in orbital character. 

In summary, we have demonstrated that in mag- 
netic systems with sufficiently large exchange splitting, 
AE 0X 3> Aejf-M, i.e. when majority (|) and minority 
spin {{) electronic states are well separated, even small 
Rashba splittings can be determined by two subse- 
quent measurements with opposite sample magnetiza- 
tions. Since one can expect an analogous behavior for 
other magnetic materials, a measurement of the Rashba 
splitting opens up a new and powerfull way to discrim- 



inate surface and interface magnetism from bulk mag- 
netism. The present discovery of a particularly large 
Rashba effect at an interface between a two-dimensional 
metal oxide and a magnetic metal should stimulate fu- 
ture research towards a potential use of such interfaces 
for spintronic devices. - Moreover, based on the present 
findings it appears to be particularly interesting to study 
the evolution of an exchange-split two-dimensional state 
into a laterally confined quantum-well state of a mag- 
netic nanostructure where the clastic reflection of the 
state |k||) is suppressed owing to the Rashba effect, since 
the reflected state | — k||) is energetically not accessible. 
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